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PROBLEMS  OF  TWO-DIMENSIONAL  ANALYSIS  OF  PARTICLE  DISTRIBUTION  IN 
AN  ATMOSPHERIC  AEROSOL 


K.  B.  Yudin 

For  studying  the  microphysical  characteristics  of  an  aerosol 
extensive  use  is  being  made  of  photoelectric  sensors,  the  operation 
of  which  is  based  on  the  measurement  of  the  light  which  is  scat- 
tered by  each  particle  at  the  moment  when  it  passes  the  illuminated 
measuring  light  cavity.  The  pulses  of  scattered  light,  which  are 
recorded  by  a photomultiplier,  are  converted  into  electrical  pulses 
and  are  fed  to  specialized  computing  and  recording  devices.  The 
dimensions  of  the  light  measuring  cavity  of  the  sensor  are  made 
quite  small  so  that  the  probability  of  the  simultaneous  entry  of 
two  particles  into  it  with  the  assigned  aerosol  concentration 
would  be  low. 

The  processing  of  the  results  of  the  measurements  Is  done 
with  the  assistance  of  the  device  of  mathematical  statistics.  In 
article  [2],  monograph  [5],  and  a number  of  other  sources  a method 
has  been  developed  for  calculating  the  physical  characteristics 
of  an  investigated  aerosol  medium  (visibility,  water  content)  and 
the  coefficients  of  the  approximating  functions  for  the  distribu- 
tion of  particles  by  size,  i.e.,  the  parameters  which  result  from 
a unidimensional  distribution.  The  contemporary  level  of  tech- 
nology makes  it  possible  to  conduct  a second  form  of  analysis  - 


1. 


to  study  the  distribution  of  time  intervals  between  the  particles 
which  are  arriving  for  measurement,  which  carry  information  on  the 
distribution  of  particles  in  space.  Two-dimensional  analysis 
makes  it  possible  to  appraise  the  degree  of  connection  of  particle 
size  with  the  distances  between  them  and  their  mutual  influence, 
and  in  particular  to  evaluate  the  fluctuations  of  mean,  mean-square 
and  mean-cubic  diameters  under  the  influence  of  spatial  distribu- 
tion. The  main  part  of  the  present  work  is  devoted  to  the  mathe- 
matical aspects  of  the  analysis  of  two-dimensional  distribution. 

The  basic  relationships  are  considered  which  describe  the 
distribution  of  a system  of  two  random  processes  relative  to  the 
problem  of  studying  the  microstructure  of  an  aerosol.  On  the 
basis  of  the  relationships  for  calculation  of  the  numerical  char- 
acteristics of  distribution  a method  is  proposed  for  evaluating 
their  fluctuations  under  the  action  of  the  spatial  distribution 
of  particles.  Since  the  analysis  of  microstructure  is  carried  out 
with  the  help  of  discrete  specialized  computers,  the  analogous 
expressions  are  given  also  in  a discrete  form.  A separate  analysis 
is  made  of  a case  in  which  the  distribution  of  particles  in  space 
is  assumed  known  and  subordinate  to  Poisson  statistics.  In  the 
concluding  part  of  the  article  the  question  of  the  selection  of 
the  Intervals  of  quantization  and  the  number  of  channels  of  the 
analyzer  In  the  study  of  particle  distribution  by  size  is  considered. 

1.  General  Relationships 

The  signal  on  the  output  of  the  photoelectric  converter  rep- 
resents the  randomly  distributed  in  time  sequence  of  electrical 
pulses  of  random  amplitude  and  can  be  described  by  the  two-dimen- 
sional integral  law  of  distribution 

F[D,  t)  Pj(D  < D,)U  < *)J. 

where  P - probability  operator,  D - particle  size,  D1  - certain 
current  value  of  particle  size,  t - random  time  Interval  between 
two  particles  arriving  in  succession  for  analysis,  x - certain 
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Then  the  function  of  distribution  can  be  written  in  the  form 

/)  T 

| | f(U,  ’)dDdt, 

^tnin  "min 

from  which  the  densities  of  distribution  of  each  of  the  processes 
are  readily  determined: 

’mux  nmtr 

f fiD,  x)dx  n /,(•*)  = j f{/),-)dO. 

’min  "min 

We  will  denote  the  conditional  densities  of  distribution  of  vari- 
ables, calculated  under  the  condition  that  the  second  variable 
accepted  an  assigned  value,  through 

?,  (DM  n h//)). 

In  order  to  avoid  errors  it  is  necessary  to  keep  in  mind  that  with 
a fixed,  for  example,  particle  size  0=0.^  the  conditional  density 
<j>2(T/Di)  characterizes  not  the  intevals  of  time  between  two  parti- 
cles of  the  same  size  which  are  arriving  in  sequence  for  analy- 
sis, but  the  intervals  of  time,  preceding  the  particles  with  the 
size  D,,  and  following  behind  the  preceding  particles  regardless 
of  their  size.  In  the  same  manner  in  the  case  of  a fixed  x=Xj 
the  conditional  density  ^(D/x^)  characterizes  the  distribution 
by  size  of  the  particles  which  are  arriving  for  analysis  after  the 
interval  Xj  after  the  preceding  particle,  and  not  the  distribution 
of  particle  size  on  assigned  intervals  of  time. 
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The  density  of  distribution  of  the  system  is  determined  by 
the  expression 

f\l),  x)— /,  (U)<t,^;/))  ./;( t)?,  (Uh). 

If  the  investigated  variables  are  independent  of  each  other  (such 
a situation,  for  example,  will  take  place  in  the  case  of  strict 
fulfillment  of  the  Poisson  law  for  the  distribution  of  particles 
in  space),  then  the  problem  is  simplified  and 

(2) 


2.  Numerical  Characteristics  of  Distribution 

In  the  course  of  the  test  the  investigators,  as  a rule,  take 
several  measurements  of  the  function  or  density  of  distribution 
of  particles.  Each  of  such  measurements  is  called  a realization, 
and  the  volume  of  aerosol  containing  a certain  number  of  particles 
which  is  taken  for  analysis  in  one  realization  - a sample.  In 
the  subsequent  account  it  is  assumed  that  the  investigated  process 
is  stationary  in  the  course  of  the  test. 

The  two-dimensional  analysis  of  particle  distribution  by  size 
and  interval  is  made  with  the  help  of  specialized  computing  devices 
- amplitude-time  analyzers.  Based  on  the  distribution  obtained  in 
one  realization  the  average,  mean  square,  and  mean  cubic  dimensions 
of  the  particles  are  calculated  using  the  formulas 


/J,=  j l)f,(D)dn,  /;,=  (/  I f)'  /,  (D)  ,/D  . 
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For  increasing  accuracy  by  means  of  excluding  errors  of  a 
statistical  nature  the  average  dimensions  are  frequently  calculated 
on  K realizations  on  the  basis  of  the  relationship 

^1=^7  2 J 

» °mln 


For  the  mean  square  and  mean  cubic  dimensions  correspondingly 


U 


D1  ft  (U)dlJ  , 


3 


n , 


|/v£ 


ir  /,  (D)  dD 


The  average  interval  between  particles  arriving  for  analysis  is 
calculated  by  the  formulas 


* = J ’/i  (T)  d~, 

’min 


j </.<<>*• 

k ’min 


The  degree  of  connection  can  be  characterized  by  the  correla- 
tion moment,  representing  the  second  mixed  central  moment  of  a 
system  of  two  random  processes  and  determined  by  the  expression 

'mu 

/?„.,=  j J (D  — 

"min 


As  is  known,  the  difference  of  the  correlation  moment  from  zero 
proves  the  dependence  of  the  processes.  Since  the  correlation 
moment,  as  this  is  evident  from  the  relationship  given,  character- 
izes not  only  the  degree  of  connection  of  the  variables,  but  also 
their  scattering,  then  for  the  evaluation  of  the  degree  of  connec- 
tion in  a "pure  form"  it  is  convenient  to  switch  to  a coefficient 
of  correlation,  equal  to 
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Ru,i 

1,,"~  8n8,  ’ 


where  and  ot  - mean  square  deviations  of  variables  D and  x, 
equal  to 


(U  — D)3f}(D)(iU 


’m«i 

f (t  - . 

*mln 


Indirectly,  apparently  the  connection  under  consideration  can  be 
judged  based  on  the  coefficient  of  correlation  between  the  average 
diameters  and  the  concentration  of  particles  in  K realizations. 

It  is  interesting  to  appraise  the  influence  of  spatial  distribution 
of  particles  on  the  average  dimensions  calculated  according  to 
formulas  (3).  For  this  purpose  it  is  possible  to  calculate  the 
dispersion  of  average  dimensions  which  is  caused  by  the  random 
distribution  of  time  Intervals  between  the  particles  which  have 
arrived  for  analysis 


3],  (T)~// 


pf  (/),  t)  an  = f ip  ji  j /(«.  t) | (in. 


i.  '■* 


Here  and  subsequently  fl,  - dispersion  operator. 

With  a calculation  that  dispersion  is  calculated  in  a point, 
equal  to  the  average  of  one  of  the  variables  which  enters  into  the 
system,  it  is  possible  to  write  [1] 


•?,,(*>-  j n*Ji\*(D)\dQ. 


'ml  n 


(4) 
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For  the  mean  square  and  the  mean  cube  dimensions  the  corresponding 
expressions  will  be 


^(•)-  ;/  r DWl\*(D)}Dcl  , 

I "mix 

»?,,(•)-  •/  f /ry/|t(/jjirfo  . 

It  should  be  noted  that  the  question  of  the  mutual  connection  of 
dimensions  of  particles  and  the  intervals  between  them,  and  moreso 
the  degree  of  connection,  has  not  been  studied  at  the  present  time. 
In  the  overwhelming  majority  of  works  on  the  study  of  the  micro- 
structure of  an  aerosol  only  a unidimensional  analysis  of  the  dis- 
tribution of  particles  by  size  is  made.  In  this  case  for  the  cal- 
culation of  the  average  diameters  expressions  (3)  are  used. 


3.  Discrete  Form 

For  measuring  the  laws  of  distribution  of  particles  with  a 
two-di.nensional  analyzer  the  magnitudes  of  frequencies  » cor- 
responding to  the  probability  of  entry  of  particles  into  each  of 
the  channels  of  the  analyzer,  are  defined  as  the  ratio  of  the  num- 
ber of  particles  n^  of  diameter  D^,  measured  in  the  given  channel, 
which  have  arrived  for  analysis  during  the  interval  of  time  Tj 
after  the  preceding  particle,  to  the  overall  number  of  all  the 
particles  recorded 


The  probability  densities  of  each  of  the  processes  can  be  found 
using  the  formulas 

7. 
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/.<a,)  = v>yi  v^. 

/ i 

In  the  case  of  unidimensional  analysis  the  magnitudes  of  frequen- 
cies are  determined  from  the  relationship 


The  average,  mean  square,  and  mean  cubic  diameters  of  one  realiza- 
tion are  calculated  using  the  formulas 


3 


The  mean  square  deviation  of  any  of  the  diameters 


,/  V 
1 r 


By  p, 


where  the  symbol  s shows  the  deviation  of  which  of  the  diameters 
is  calculated.  For  processing  the  results  of  the  tests  on  K 

realizations  it  is  convenient  to  introduce  the  following  relation- 
ships : 


fh 


/ -i 


* 


/-i 


* 


then 


B = p,  . Dt  = 1 / V pt  , /;,  = ./  V D,pi  , 

I • / • « 


°/»,  = 1/ 


/ 1 V,r, 


For  the  intervals  of  time  between  particles 

tr:  •'/'/•  3,  = 1 (*/ - ")"/'/  • 

/ 

The  correlation  moment  is  calculated  using  the  relationship 

/v’ot  — ^ 2^  ^ i — (Ty  Ptj  • 

/ » 

The  correlation  moment  between  average  diameter  of  particles  and 
their  concentration  in  K realizations  can  be  found  from  the 
relationship 

A * 

^T>n  ~ ^ (^1  fi)  (^  / ~ W). 

(-1  /-i 

Here  the  doubled  sum  on  1 indicates  that  all  the  products  in  the 
case  of  all  values  of  1 for  each  of  the  factors  should  be  summed. 
The  mean  square  deviations  are  found  using  the  formulas 


l ~ I 


The  dispersion  of  average  diameter  under  the  influence  of  spatial 
distribution  of  particles  is  equal  to 


V Dt  p,  </).-) 


V/;j,V  (D,  t)|  = V0?/riA(*>l- 

i l 


For  the  dispersions  of  the  mean  square  and  mean  cubic  diamters 
the  following  expressions  are  valid 


L 
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In  the  expression  for  average  diameters,  calculated  on  K realiza- 
tions, instead  of  A[pi  (^ ) 3 we  introduce  AEp^t)]. 

The  coefficients  of  variations  (variability)  of  average  diam- 
eters in  all  cases  are  calculated  as  the  ratio  of  the  mean  square 
deviation  to  the  corresponding  average 


4.  Distribution  According  to  Poisson  Law 

The  evaluation  of  a on  the  Dasis  of  expression  (4)  can  be 
done  with  unidimensional  analysis  based  on  the  a priori  known 
form  of  distribution  of  particles  in  space.  According  to  concepts 
which  exist  at  the  present  time  this  distribution  is  subordinate 
to  Poisson  statistics.  In  this  case  the  probability  of  finding 
n^  particles  of  diameter  in  a unit  measured  volume  is  given 
by  the  expression 


Pi  = 


"1 

«il 


This  distribution  has  mathematical  expectation  and  dispersion, 
equal  to  n^.  Since  the  concentration  of  particles  of  different 
diameters  is  different,  i.e.  fti=f(D),  then  it  is  evident  that  the 
dispersion  of  intervals  between  particles  of  different  diameters 
is  different  and  is  proportional  to  their  concentration,  i.e. 

//|t  (0)1 -»(«). 

The  last  expression  shows  that  the  dispersion  of  average  diameters 
can  be  calculated  only  on  the  basis  of  several  realizations,  and 
their  number  should  be  sufficient  enough  so  that  the  evaluation 
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of  n(D)  would  be  significant. 

Having  substituted  the  value  A[x(D)]  Into  the  expression, 
for  example,  of  dispersion  of  average  diameter  under  the  influence 
of  the  spatial  distribution  of  particles,  in  the  case  of  distri- 
bution by  Poisson  law  we  obtain 


^)"T„  = ~ j l>'n{D)dD. 


(5) 


nm\n 


In  the  same  manner  the  expression  for  mean  square  and  mean  cubic 
diameters  is  obtained. 

In  a discrete  form,  taking  into  account  that  /»,  = —'  , we 

X 


obtain 


and  for  K realizations 


A,J  A'M  AW*  A'/V 


Dispersion  and  the  coefficient  of  variations  of  average  di- 
ameter in  K realizations  in  the  case  of  Poisson  distribution  of 
particles  in  space  will  be  equal  to 


ft 


(6) 


1 
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These  parameters  can  also  be  calculated  using  experimental  data 
and  the  following  relationships: 


T 


i rs 


(5„ 


0,)J.  r/ie  Dt  =-jJ-S  7J,> 

k 


And  if  it  turns  out  that  3r».('t) ny.r  ^ S/T , then  the  fluctuations 

of  D-^  can  be  explained  by  fluctuations  in  the  distribution  of 
particles  in  space.  When  ’u  n,K^  ’n,™,!  and  correspondingly 

ii v«  < '*7>, oMMt  , it  follows  to  assume  that  the  fluctuation  of 
average  diameter  is  connected  with  fluctuations  in  the  distribu- 
tion of  particles  by  size.  What  was  said  is  also  valid  for  the 
mean  square  and  mean  cubic  diameters  and  the  characteristics  of 
an  aerosol  which  are  calculated  on  their  basis. 

Let  us  consider,  for  example,  the  fluctuation  in  water  content 
of  an  aerosol  with  an  unchanged  distribution  of  particles  by  size. 
The  water  content  of  the  aerosol  is  calculated  using  the  formula 


* * 


The  dispersion  of  water  content  in  K realizations,  caused  by  the 
fluctuation  in  the  number  of  particles  due  to  their  random  dis- 
tribution in  space,  will  be  determined  by  the  expression 


In  the  case  of  distribution  of  particles  in  space  according  to 
Poisson  law  this  expression  acquires  the  form 


4 llyw 


/r,v  y 7 


, 2^* 


KX 
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i.e.,  dispersion  is  proportional  to  N.  The  corresponding  coeffic- 
ient of  variatior 


* 


I • 


5.  Intervals  of  Quantization  and  Number  of  Channels  for  the 
Analyzer 


An  important  methodical  problem  in  the  study  of  the  distri- 
bution of  particles  by  size  is  the  rational  selection  of  the  inter- 
vals of  quantization  and  the  number  of  channels  for  the  threshold 
analyzer.  Below  we  will  consider  one  of  the  possible  methods  of 
evaluating  these  parameters  which  is  based  on  knowledge  of  the  law 
of  measured  distribution  and  permissible  error.  For  this  purpose 
the  density  of  probability  pi=f1(D1)  is  approximated  by  a poly- 
nomial of  the  first  degree  for  each  interval  Di+1-Di  (piecewise 
linear  approximation).  The  error  of  such  approximation  is  determ- 
ined by  the  remainder  term  of  the  interpolation  formula.  For  the 
Newton  formula  the  maximum  magnitude  of  the  remainder  term  is 
equal  to  [4] 


I J 


fin,)  /Dm  — o,  y 
' 2!  \ 2 ) 


Assuming  the  absolute  error  Ap.^  nax=max  I J(Di ) I > for  determination 
of  the  sector  of  partitioning  AD*(Di+1-D)  in  [3,  6]  the  following 
expression  is  obtained 


M)~ 


13. 


where  |f”(D)|max  - maximum  absolute  value  of  the  second  deriva- 
tive in  point  D of  the  interval  D..,-D.  or  for  relative  error 

l+J.  1 


f\fJ)  ' 


M)~  I, ^8« 


rrm 


The  required  number  of  channels  for  the  analyzer  with  a uniform 
partitioning  by  particle  size  will  be  equal  to 

Dm»x  " /I min 
m ' . 

w 


For  example,  for  gamma-distribution  with  a whole  and  positive 


/> 

[)'- 2 { I \2 

no)  = -5^rr  — T°)  -«  • 


nu) w 

“F tT 
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With  the  parameters  of  gamma-distribution  a=8  and  0=1.2  which  are 
encountered  frequently  in  clouds  the  second  derivative  Js  maximum 
in  a point  which  is  close  to  4.5  pm.  If  the  range  of  measured 
dimensions  of  particles  lies  within  the  limits  from  2 to  30  pm, 
then  with  a relative  error  e*0.05  we  have 
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= 2.3.  A0  = O.%.  m 


quantization  in  all  the  remaining  inter- 
vals will  be  less  than  assigned.  If  with  the  assigned  accuracy 
it  is  required  to  determine  the  distribution  in  the  area  of  aver- 
age size  D=8 (a+1 )=10 . 8 pm,  then 

pi* 
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_ — — 16  7 A/l  = 2,6  mum  , m — 12, 
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but  the  error  in  the  area  of  sizes  of  4.5  pm  increases  to  0.36. 

With  nonuniform  quantization  the  magnitudes  of  the  intervals 
AD  are  not  equal  to  each  other,  but  change  according  to  a selected 
law.  It  is  convenient  to  consider  the  advantages  of  nonuniform 
quantization  in  an  example  of  quantization  according  to  logarith- 
mic law.  It  Is  known  that  gamma-distribution  for  the  description 
of  distribution  by  size  of  cloud  particles  is  introduced  as  the 
approximating  function  for  logarithmically  normal  distribution. 
Having  designated  u=lnD,  we  have 

* = L / 
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The  thresholds  of  the  analyzer  channels,  expressed  In  particle 
dimensions  (pm),  and  not  their  logarithms,  should  be  set  on  the 
levels  2.0,  2.9,  4.5,  6.6,  9-9,  14.8,  22.0,  32.9. 


Thus  due  to  the  expansion  of  channels  in  the  area  of  large 
particles  the  required  number  of  channels  was  reduced  to  7,  and 
the  error  from  quantization  was  not  increased. 
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